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ABSTRACT 


We  show  how,  as  a  result  of  the  strong  heating  produced  at  chromospheric 
levels  during  a  solar  flare  burst,  the  local  gas  pressure  can  transiently 
attain  very  large  values  in  certain  regions.  The  effectiveness  of  the  sur¬ 
rounding  magnetic  field  at  confining  this  high  pressure  plasma  is  therefore 
reduced  and  the  flaring  loop  becomes  free  to  expand  laterally.  In  so  doing  it 
may  drive  magnetic  field  lines  into  neighboring,  non-flaring,  loops  in  the 
same  active  region,  causing  magnetic  reconnection  to  take  place  and  triggering 
another  flare  burst.  The  features  of  this  interacting  loop  model  are  found  to 
be  in  good  agreement  with  the  energetics  and  time  structure  of  flare- 
associated  solar  hard  X-ray  bursts. 


Subject  headings:  hydromagnetics  -  plasmas  -  Sun  :  flares  -  Sun  :  X-rays 


I.  INTRODUCTION 


Solar  hard  X-ray  bursts  (photon  energy  e  i,  10  keV)  frequently  exhibit  de¬ 
tailed  time  structure  on  timescales  of  a  few  seconds  (see,  e.g.,  Hoyng,  Brown, 
and  van  Beek  1976;  Dennis,  Frost,  and  Orwig  1981).  This  has  led  to  the  con¬ 
cept  of  the  "Elementary  Flare  Burst"  ( EFB  -  van  Beek,  de  Feiter,  and  de  Jager 
1974;  de  Jager  and  de  Jonge  1978),  whereby  it  is  supposed  that  hard  X-ray 
bursts  with  a  complicated  time  structure  involving  several  "spikes"  are  in 
fact  composed  of  a  number  of  discrete  EFB's,  each  with  a  simple  "rise-fall" 
time  structure. 

de  Jager  and  de  Jonge  (1978)  were  the  first  to  point  out  that  there  ap¬ 
pears  to  be  a  significant  difference  (as  measured  by  their  FWHMs)  amongst 
EFB's  in  a  particular  hard  X-ray  event,  suggesting  that  different  EFB's  cor¬ 
respond  to  the  successive  activation  of  a  number  of  different  sources,  rather 
than  the  continual  reactivation  of  a  single  source.  Karpen,  Crannell,  and 
Frost  (1979)  also  studied  the  structure  of  hard  X-ray  bursts  exhibiting  clear 
EFB  structure  (in  their  nomenclature,  "multiply  impulsive"  events),  together 
with  microwave  data  on  the  same  bursts,  and  showed  that  the  source  parameters 
(density,  temperature,  magnetic  field,  etc.)  for  each  EFB  were  markedly  dif¬ 
ferent.  This  similarly  led  them  to  the  conclusion  that  different  EFB's  origi¬ 
nate  in  different  regions  of  the  flare.  Finally,  Kane,  Pick,  and  Raoult 
(1980)  studied  the  spatial  structure  of  type  III  radio  bursts  observed  to  be 
synchronous  with  hard  X-ray  EFB's  and  deduced  thac  each  of  the  hard  X-ray 
bursts  in  the  event  in  question  (on  2  September  1978)  originated  in  one  n*  two 
spatially  distinct  regions,  separated  by  some  7  minutes  of  arc,  or  3  x  1010  cr.. 


There  thus  appears  to  be  a  significant  amount  of  evidence  for  different  EFB's 
corresponding  to  excitation  of  different  parts  of  the  flare  region. 

Kane,  Pick  and  Raoult  (1980),  in  attempting  to  explain  their  inferred 
spatial  separation  of  the  individual  hard  X-ray  bursts,  rule  out  chance  coin¬ 
cidence  of  the  two  bursts  occurring  independently  as  being  statistically  im¬ 
probable1,  and  point  out  that  any  signal  from  one  region  to  the  other  would 

1  Kane,  Pick,  and  Raoult  (1980)  quote  an  occurrence  rate  of  "  10  day-1  foe  the 
class  of  burst  appropriate  to  their  observations,  and  so  deduce  the  probabil¬ 
ity  of  two  bursts  occurring  within  a  five  second  time  interval  (the  temporal 
separation  of  the  bursts)  to  be  <  10“7.  This  is  in  fact  the  probability, 
based  on  Poisson  statistics,  of  two  bursts  occurring  in  a  given  5  second 
time  interval.  However,  since  the  time  of  the  first  burst  is  arbitrary,  one 
should  actually  calculate  the  probability  of  a  single  burst  occurring  in  a 
given  5  second  interval;  this  evaluates,  using  the  same  occurrence  rate 
quoted  by  Kane,  Pick,  and  Raoult,  to  be  *  5  x  10"1,  or  about  1  in  2000. 
Although  this  is  still  a  small  number,  it  is  conceivable,  considering  the 
large  number  of  bursts  observed,  that  chance  coincidence  may  in  fact  have 
been  the  responsible  agent  in  the  event  studied  by  Kane,  Pick,  and  Raoult. 

have  to  travel  much  faster  than  any  reasonable  Alfven  speed,  thus  ruling  out  a 
sympathetic  trigger.  They  thus  interpret  their  observations  as  implying  acti¬ 
vation  of  two  different  source  regions  by  the  same  external  disturbance,  i.e. 
an  energetic  electron  beam  in  a  "thin  target"  (Datlowe  and  Lin  1973)  scenario. 
However,  in  the  event  studied  by  Kane,  Pick,  and  Raoult  (1980),  the  two  flar¬ 
ing  regions  are  separated  by  a  distance  comparable  to  a  solar  radius,  and 
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therefore  clearly  belong  to  different  active  regions;  for  two  hard  X-ray 
sources  located  in  the  same  active  region,  a  sympathetic  trigger  cannot  be 
ruled  out  by  such  disturbance  velocity  considerations  (see  Karpen,  Crannell, 
and  Frost  1979).  In  this  paper  we  propose  a  means  of  producing  just  such  a 
trigger. 

To  date,  all  models  of  the  response  of  the  solar  atmosphere  to  a  flare 
energy  input  have  been  one-dimensional  in  nature.  This  simplifying  character¬ 
istic  is  usually  justified  by  noting  that  for  typical  parameters  in  the  flar¬ 
ing  corona  and  chromosphere  the  quantity  8,  defined  as  the  ratio  of  gas  to 
magnetic  pressures,  is  much  less  than  unity,  so  that  the  plasma  is  constrained 
to  move  along  the  magnetic  field  lines.  (The  assumption  of 8  <  1  breaks  down 
near  the  photosphere  [number  density  n  >  10  cm  ] ,but  such  deep  layers  of 
the  atmosphere  remain  relatively  undisturbed  during  the  flare  [Machado, 

Emslxe,  and  Brown  1978].)  However,  the  relation  8  «  1  follows  from  considera¬ 
tion  of  steady-state  conditions;  as  we  shall  show  below  (§11)  the  fact  that 
the  timescales  for  impulsive  heating,  and  for  hydrodynamic  relaxation,  of  the 
solar  chromosphere  are  markedly  different  can  result  in  the  formation  of  re¬ 
gions  of  transiently  very  large  gas  pressure,  in  which  8  is  comparable  with 
unity.  While  to  explore  fully  the  resulting  three-dimensional  time-dependent 
response  of  the  solar  atmosphere  to  the  flare  energy  input  is  beyond  the  scope 
and  purpose  of  the  present  paper,  we  sh*ll  discuss  semi -quantitatively  the  ef¬ 
fect  of  creating  such  a  transient  high-8  regime.  In  particular,  we  shall  dem¬ 
onstrate  (§111)  how  this  high  gas  pressure  region  causes  the  energized  flare 
loop  to  expand  laterally,  dragging  the  (frozen-in)  magnetic  field  lines  with 
the  plasma.  A  magnetohydrodynamic  disturbance  is  thus  established,  which  can 
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subsequently  interact  with  the  field  lines  of  a  neighboring  quasi-static  loop 
in  the  same  solar  active  region,  causing  magnetic  reconnection  and  further  re¬ 
lease  of  energy,  thus  triggering  the  next  EFB.  This  process  can  be  continued 
amongst  the  various  loops  in  the  preflare  active  region  complex,  producing  the 
"multiple-spike"  X-ray  flux  versus  time  profile  characteristic  of  so  many 
events  (Hoyng,  Brown,  and  van  Beek  1976;  Karpen,  Crannell,  and  Frost  1979). 

The  energetics  and  timescales  of  this  process  are  explored  in  §111,  and  it  is 
shown  that  both  are  compatible  with  their  counterparts  inferred  from  hard  X- 
ray  data.  In  §IV  we  state  our  conclusions  and  suggest  areas  for  more  detailed 
research. 


II.  FORMATION  OF  REGIONS  OF  TRANSIENTLY  HIGH  GAS  PRESSURE 
DURING  THE  IMPULSIVE  HEATING  PHASE  OF  A  FLARE  BURST 


There  is  some  controversy  as  to  the  exact  mechanism  responsible  for  heat¬ 
ing  the  solar  chromosphere  during  flares,  candidates  including  non-thermal 
electron  bombardment  (Brown  1973;  Lin  and  Hudson  1976;  Emslie  1980),  soft  X- 
ray  irradiation  (Somov  1975;  Henoux  and  Nakagawa  1977,  1978;  Machado  1978; 
Hyder  1981),  and  dissipation  of  hydrodynamic  (Craig  and  McClymont  1976)  and 
thermal  (Smith  and  Harmony  1981)  shocks.  For  the  purposes  of  the  present  dis¬ 
cussion  we  shall  concentrate  on  the  first  of  these,  for  two  reasons.  First, 
our  discussion  will  relate  to  events  with  a  strong  hard  X-ray  signature;  to 
date  all  models  (both  "thermal"  and  "non-thermal";  see  Emslie  and  Rust  1979) 


of  the  hard  X-ray  production  mechanism  involve  a  substantial  amount  of  the 
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flare  energy  to  be  in  the  form  of  a  beam  of  precipitating  high  energy  elec¬ 
trons  (see,  e.g..  Brown  1971;  Melrose  and  Brown  1976;  Brown,  Melrose,  and 
Spicer  1979;  Vlahos  and  Papadopoulos  1979;  Emslie  and  Vlahos  1980;  Smith  and 
Brown  1980).  Second,  and  more  importantly  for  our  present  study,  this  mecha¬ 
nism  of  chromospheric  heating  has  been  extensively  studied  in  the  literature; 
in  particular,  calculations  of  the  detailed  time-dependent  response  of  the 
solar  atmosphere  to  such  a  beam  of  non-thermal  electrons  have  been  carried  out 
by  a  number  of  authors  (Kostyuk  and  Pikel'ner  1975;  Sermulina  et  al.  1980; 
Somov,  Syrovatskii,  and  Spektor  1981). 

Early  work  on  the  response  of  the  solar  chromosphere  to  a  hypothetical 
flare  energy  input  was,  for  reasons  of  simplicity  and  tractability ,  confined 
to  a  discussion  of  steady-state  model  atmospheres  (e.g.,  Brown  1973;  Henoux 
and  Nakagawa  1977,  1978;  Brown  Canfield,  and  Robertson  1978;  LaBonte  1978). 
Kostyuk  and  Pikel'ner  (1975)  attempted  a  time-dependent  calculation  of  the 
hydrodynamic  response  of  the  solar  atmosphere  to  a  sustained  (»  100  s)  injec¬ 
tion  of  non-thermal  electrons.  However,  they  did  not  adequately  treat  the 
problem  of  radiative  instability  of  the  heated  plasma;  further,  their  pub¬ 
lished  results  lack  important  information,  such  as  density  profiles  (see 
Canfield  et  al.  1980).  Recently  more  detailed  work  on  the  time-dependent  re¬ 
sponse  of  the  atmosphere  to  an  injected  beam  of  non-thermal  electrons  has  been 
carried  out  (Sermulina  et  al.  1980;  Somov,  Syrovatskii,  and  Spektor  1981). 

Due  to  numerical  difficulties,  these  analyses  are  restricted  to  only  short 
<  10  second)  bursts  of  heating,  and  subsequent  relaxation  of  the  energized  at¬ 
mosphere  to  its  pre-flare  state.  This  timescale  is,  however,  quite  acceptable 
for  studying  events  composed  of  discrete  EFB's  (de  Jager  and  de  Jonge  1978; 
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Figure 


Karpen,  Crannell,  and  Frost  1979).  Further,  the  results  clearly  show  that  in 
all  cases  the  electron  temperature  of  the  pre flare  chromospheric  material  is 
heated  to  coronal  temperatures  in  times  on  the  order  of  one  second,  while  the 
density  of  this  material  stays  relatively  constant,  due  to  the  much  longer  hy¬ 
drodynamic  time  scale.  Knowledge  of  the  detailed  behavior  of  the  atmosphere 
long  after  this  impulsive  rise  in  electron  temperature  is  not  necessary  for 
the  subsequent  discussion. 

Figure  1  (from  Somov,  Syrovatskii,  and  Spektor  1981)  shows  the  density 

n  (cm~^)  and  electron  temperature  Tg  (K)  as  a  function  of  column  density  N 

_  2 

(cm  )  after  5  seconds  of  heating  by  a  beam  of  non-thermal  electrons  with  an 
injected  energy  flux  spectrum  (electrons  cm  ‘‘s  per  unit  energy)  in  the  form 
of  a  truncated  power  law: 


where  6=3,  =  10  keV,  and  5  (the  total  injected  energy  flux)  =  1011  ergs 

cm  s  The  behaviors  of  n  and  Tg  in  the  (hydrostatic)  preflare  atmosphere 
are  shown  dashed.  Note  that  in  this  preflare  state  Te  is  constant  (  =  6700  K) 
and  n,  consequently,  has  an  exponential  variation  with  height  (so  that  n  <r  N). 
This  initial  state  is  admittedly  somewhat  unphysical;  however,  it  turns  out 
(Somov,  Syrovatskii,  and  Spektor  1981)  that  the  detailed  structure  of  the 
preflare  atmosphere  at  low  column  densitites  (great  heights)  is  not  an  impor¬ 
tant  consideration.  This  is  because  at  such  low  column  densities  the  flare 
energy  input  is  much  larger  than  any  preflare  energy  source  or  sink  (see 
Emslie  1980),  so  that  the  plasma  is  explosively  heated,  resulting  in  a  final 
structure  which  is  essentially  independent  of  the  assumed  initial 


conditions.  At  greater  column  densities  (below  the  region  of  explosive 
heating  —  see  Brown  1973),  the  preflare  and  flare  energy  inputs  are  more 
comparable,  and  due  consideration  must  be  given  to  the  initial  conditions  of 
the  atmosphere.  Somov,  Syrovatskii,  and  Spektor's  (1981)  assumed  value  of  T0 
=  6700  K  reflects  this  point,  since  this  is  a  typical  temperature  for  the 
upper  chromosphere  (see  Machado  et  al.  1980).  Clearly  a  more  detailed 
description  of  Tg(N)  for  large  N  would  be  desirable  in  order  to  more  accurate¬ 
ly  assess  the  structure  of  the  flaring  chromosphere  at  these  levels.  However, 
this  paper  will  be  concerned  only  with  the  structure  of  the  atmosphere  in  the 
explosively  heated  region,  and  so  we  shall  not  dwell  on  this  matter  further. 
Also  shown  in  Figure  1  is  the  flaring  pressure 

P  =  nTe  cm-3  K  (2) 

and  the  equipartition  magnetic  field  strength 

Beq  =  <8wkP)1/2  G  '  (3) 

where  k  is  Boltzmann's  constant.  In  writing  eq.  (2),  we  have  neglected  the 

contribution  from  the  ions,  since  T.  «T  (Somov,  Syrovatskii,  and  Spektor 

l  e 

1981).  It  is  clear  from  the  figure  that  a  region  of  greatly  enhanced  gas 

pressure  is  formed  in  the  range  lO1^  <  N  <  4  x  10  ^  cm-2.  For  101g  <  N  <  2  x 

10  cm  this  is  the  result  of  a  sudden  increase  in  temperature  over  preflare 

1 9 

values  (the  density  remaining  essentially  unchanged),  while  for  2  x  10 
-2  19  —2 

cm  <  N  <  4  *  10  cm  the  excess  pressure  is  attributable  to  a  density  en¬ 
hancement,  brought  upon  by  compression  due  to  the  overlying  high  pressure 
material.  This  region  of  enhanced  density  is  very  thin,  being  only  AN/n  =  2 
km  thick,  and  is  a  characteristic  feature  of  the  chromospheric  side  of  a 
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(transient)  flare-induced  transition  region.  When  heated  it  expands  at  con¬ 
stant  pressure  to  form  a  region  with  the  original  preflare  density  and  a  new 
hotter,  temperature,  and  a  new  region  of  enhanced  density  is  created  somewhat 
deeper  in  the  atmosphere.  As  the  heating  proceeds,  this  layer  (and  so  the 
transition  region)  "burns"  its  way  deeper  and  deeper  into  the  atmosphere. 

This  process  stops  when  the  heat  being  deposited  on  the  chromospheric  side  of 
this  region  has  been  reduced,  by  attenuation  of  the  electron  flux  by  the 
overlying  material,  to  such  an  extent  that  it  can  no  longer  overcome  the 
ability  of  the  plasma  to  radiate  the  energy  away.  The  plasma  is  now  unable  to 
enter  the  temperature  regime  of  radiative  instability  (see  Field  1965;  Cox  and 
Tucker  1969)  and  so  heat  up  to  coronal  temperatures  ( *  107  K).  The  atmosphere 
thus  remains  approximately  in  this  state  (ignoring  heat  redistribution  due  to 
thermal  conduction  and  radiative  transfer)  until  the  cessation  of  the  electron 
energy  input,  at  which  point  it  starts  to  relax  back,  through  both  thermal 
conduction  and  radiation,  to  its  preflare  state. 

Although  the  results  of  Somov,  Syrovatskii,  and  Spektor  (1981)  are  rather 
specific,  we  may  use  the  physical  insight  gained  in  the  above  discussion  to 
generalize  their  results  to  other  cases.  The  level  of  the  transition 

region  will  be  determined  by  the  total  amount  of  energy  available  for  heating 
material  to  coronal  temperatures,  while  the  magnitude  of  the  gas  pressure  in 
the  transient  regime  at  this  newly-formed  transition  region  will  be  approxi¬ 
mately  given  by 


P 


P 

o 


(4) 
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where  PQ  is  the  initial  pressure  at  N  =  NTR,  Tc  and  T  ^  are  the  coronal  and 
chromospheric  electron  temperatures  (”  107K  and  10^  K)  respectively,  and  the 
factor  2  in  the  denominator  follows  from  the  fact  that  Tg  >  during  the  im¬ 
pulsive  phase  of  the  flare,  while  Tg  a  in  the  preflare  chromosphere. 

We  now  proceed  to  estimate  HpR  from  simple  theoretical  considerations. 

The  heating  rate  (ergs  cm-^s-^)  due  to  Coulomb  collisions  of  a  beam  of  non- 
thermal  electrons,  injected  vertically  with  energy  spectrum  Fo(Eq)  is  (Brown 
1973;  Emslie  1978 — we  ignore  the  effect  of  reverse  current  ohmic  losses  in  the 
impulsive  phase  [see  Emslie  1980,  1981]) 

F  (E  )dE 

_ £_£ _ 2 _  fS1 

2  2/3  '  S 

E  (1-3KN/E  )  7 
o  o 

4 

where  K  =  2ite  A  (e  =  electronic  charge  [e.s.u.],  A  =  Coulomb  logarithm),  and 

E*  =  max  (E  ,  [3KN] 1/2)  ,  (6) 

E.j  being  the  low  energy  cutoff  to  the  injected  beam  energy  spectrum — see  eq. 
(1).  The  total  heating  rate  (ergs  cm  s  )  down  to  column  depth  N  is  thus 

N  N  I  (N’)dN’ 

JB(N)  =  /  Ib(N')  dN’  =  /  - n -  '  (7) 

0  0 

or,  using  eqs.  (1)  and  (5), 


I  (N)  =  Kn  / 
B  E- 


J0(N)  =  K  (6-2)  5  E, 


6-2 


N 


/  / 


dE  dN’ 
o 


N ' =0  E  =E *  E  6  +  1 ( 1-3KN '/E2 )2/3 
o  o  o 


(8) 
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2 

For  N  <  N1  =  Ej/3K,  we  have  E*  =  E^;  reversing  the  order  of  integration  then 
yields 


VN)I 


=  y  [i  -  ( 6-2 )e 


6-2 


N<N. 


E1_^( 1- 
o 


3KN/E2 ) 1/3dE  1 

O  O  ' 


(9) 


Setting  EQ  =  E^  sec  0  reduces  JB(N)  to  a  form  convenient  for  computation: 


J  (N)|  =  f  [l 

‘  N<N 


*/2  . 

(6-2)  /  cos  0  (1 
o 


2m1/3  • 
cos  0)  sin 


0  d0] 


(10) 


where 


3KN  __  N_ 

„2  *  N . 


(ID 


For  N  E*  =  ( 3KN )^2  :  thus,  in  evaluating  Jg(N)  for  N>N1  we  must  split  the 

integral  over  N'  into  two  parts — one  from  N'=0toN’=N1  and  the  other  from 
N'  =  N -j  to  N'  =  N.  Hie  first  part,  Jg-j,  is  given  by  equation  (10)  with  v=  1, 


TB<!  -  T>1  - 


(12) 


where  B  is  the  Beta  function  (see  Abramowitz  and  Stegun  1965,  p.  258).  Hie 
second  term, 


JB2(N)  =  (6-2)JeJ‘2  /  / 


N 


dE  dN' 
o 


N'=N  E  =(3KN')1/2  E5+1  (l-3KN’/E2  )2/3 
1  o  o  o 


(13) 
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can  be  evaluated  in  a  straightforward  manner  to  give  (see  Emslie,  Brown,  and 
Donnelly  1978) 


jb2(N)  -  i B  21 


(14) 


Finally,  adding  equations  (12)  and  (14)  gives 


i-  i 

Jg(N)|  =  5[l-  j  (n/N^  2]  ; 


(15) 


N>N 


note  that  this  result  can  also  be  obtained  by  using  the  identity 


«  I  (N* )dN ' 

?  .  J  (»)  =  Jn(N)  ♦  /  - - - 


(16) 


where  the  second  term  is  evaluated  in  a  straightforward  manner  with  E*  = 
(  3KN  '  )^2  . 


Equations  (10)  and  (15)  are  the  results  needed  to  continue  our  analysis. 
In  the  case  6  =  4  (a  not  unreasonable  value — see  Hoyng,  Brown,  and  van  Beek 
1976)  they  reduce  to  the  simple  expressions 


Jb(N)  = 


3N 


3N 


;  N  >N 


(17) 


We  shall  use  these  expressions  hereafter.  The  effect  of  varying  5  will  be 
considered  briefly  below. 
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Now  consider  a  beam  injected  for  a  time  r,  which  deposits  its  energy 
above  column  depth  HpR  and  heats  this  material  from  chromospheric  to  coronal 
temperatures •  Energy  balance  in  the  region  N  <  N^,R  thus  dictates  that 


W  T 


Nk(T  -T  > 
TR  c  ch 


N^ltT  , 
TR  c 


(18) 


where  we  have  ignored  energy  losses  due  to  thermal  conduction  and  radiation  (ci.  above); 

this  neglect  will  cause  the  NTR  calculated  from  eq.  (18)  to  be  an  upper  limit. 

19  -2 

For  E1  =  10  keV  (see  Somov,  Syrovatskii,  and  Spektor  1981),  N1  =  10  cm 
(Emslie  1978);  thus,  by  appealing  to  Somov,  Syrovatskii,  and  Spektor’s  (1981) 
results,  we  may  take  NTR  >  N1  (we  shall  in  fact  verify  this  assumption  a^ 
posteriori ) .  Using  eq.  (17)  ineq.  (18)  then  yields 


(19) 


where 


x  = 


4N 


TR 


3N  < 


;  C  = 


N.kT 
1  c 

Jr 


Equation  (19)  solves  to  give 


(20) 


x  =  [l  +  d-3C)1/2  ]  . 


(21) 


where  the  positive  sign  in  the  square  brackets  has  been  chosen,  so  that  as 
J ♦*  (C+0),  N  >  N^.  In  terms  of  physical  variables,  eq.  (21)  is 


MMHI 


i 
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V  = 


JLL 

2kT 


{l  +  [  1 -3N 1  kT^/JT  j'/^}  . 


(22) 


Now,  by  eq. 
by 


(4),  the  transient  (i.e.  flaring)  pressure  at  this  level  is  given 


P 


Nm_m„g  T 
TR  H^®  c 


2kT 


ch 


(23) 


where  mR  is  the  hydrogen  mass  and  ga  the  solar  gravity.  Substituting  for  NTR 
from  eq.  (22)  and  inserting  numerical  values  gives 

2 

<7  41 . 4E  ,T 

P  »  6  x  1011  p-  {l+[l - ~]  /2  }  •  (24) 

ch 


We  may  now  compare  this  expression  with  the  results  of  Somov, 

Syrovatskii,  and  Spektor  (1981);  see  Figure  1.  Ignoring  the  discrepancy  in 
the  value  of  6  used  (for  other  6  in  the  range  3-6  the  integrals  in  eqs.  [10] 
and  [13]  may  be  evaluated  numerically  [or  analytically]  to  yield  results  which 
differ  from  those  for  6  =  4  by  only  some  10%),  we  find  from  eq.  (22),  with  J  * 
1011  ergs  cm~2s-1,  E1  =  in  keV,  t  =  5s,  Tch  =  104  K,  and  Tc  =  107  K,  that 


N 


TR 


,  „  ,„20  -2 

1.8  x  10  cm 


(25) 


(note  that  this  is  indeed  greater  than  N1(  as  was  earlier  assumed).  Further, 
eq.  (24)  gives 


P  *  6  x  10  ^cm 


(26) 


corresponding  to  (eq.  [3]) 


B  *  450  G  .  (27) 

eq 

Although  these  values  are  somewhat  larger  than  the  values  given  by  Somov, 
Syrovatskii,  and  Spektor  (1981)  (see  Figure  1),  we  note  that  the  value  of  NTR 
in  eq-  (25)  is  more  compatible  with  empirical  estimates  of  this  quantity 
(Machado  et  al.  1980),  possibly  indicating  deficiencies  in  the  calculations  of 
Somov,  Syrovatskii,  and  Spektor  (Emslie,  Brown,  and  Machado  1981).  Even  allow¬ 
ing  for  the  fact  that  our  estimates  of  are  necessarily  upper  limits,  due 

to  the  oversimplistic  energy  balance  equation  (18),  it  nevertheless  appears, 

therefore,  that  the  gas  pressure  P  can  easily  attain  transient  values  of  order 
19  -3 

10  '  cm  K  during  the  impulsive  phase  of  flares,  correponding  to  an  equipar- 
tition  magnetic  field  strength  B_  "  200  G.  (For  example,  increasing  t  to  10s 

eM 

in  the  calculations  of  Somov,  Syrovatskii,  and  Spektor  [1981]  should  increase 

V3 

Bec(  by  a  factor  *  2  *■  over  the  values  in  Figure  1.) 

These  values  of  B  are  not  small  compared  to  the  ambient  field  strength 

eq 

in  the  chromosphere.  Thus  it  appears  that  the  assumption  of  a  plasma  con¬ 
strained  to  move  along  preexisting  field  lines  is  not  adequate  in  describing 
the  impulsive  phase  of  a  solar  flare  burst,  since  there  is  no  longer  an  over¬ 
whelming  force  (due  to  magnetic  pressure)  restricting  motion  perpendicular  to 
chese  field  lines.  In  the  next  section  we  shall  discuss  the  implications  of 
this  situation  for  the  evolution  of  the  flare. 
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III.  THE  INTERACTION  OF  NEIGHBORING  LOOPS 


2 

Consider  now  the  situation  depicted  in  the  left  sketch  of  Figure  2,  which 
shows  two  typical  loops  in  an  active  region  loop  complex  with  a  simple  bipolar 
magnetic  field  topology.  (For  clarity,  the  field  lines  between  and  around 
loops  Lj  and  Lj  are  not  shown.)  At  some  instant  loop  L^  becomes  unstable  (due 
to,  for  example,  the  formation  of  multiply  connected  magnetic  islands  in  the 
loop  [Spicer  1976,  1977]  or  the  emergence  of  new  flux  at  its  base  [Heyvaerts, 
Priest,  and  Rist  1977] )  and  flares,  as  indicated  by  the  star  at  the  loop  apex. 
This  causes  a  burst  of  hard  X-ray  emission  to  occur.  Energy  is  transported 
along  the  loop  and  deposited  at  the  dense  footpoints,  causing  a  transient 
strong  pressure  buildup  (§11).  Since  the  solar  atmosphere  is  highly  conduct¬ 
ing,  the  field  and  plasma  are  "frozen"  together,  and  so  the  large  pressure 
gradient  between  the  inside  and  outside  of  L.  drives  a  lateral  expansion  of 
the  field  lines  defining  loop  .  (Of  course,  there  are  in  addition  motions 
along  loop  L.,  driven  by  longitudinal  pressure  gradients  [see  Craig  and 
McClymont  1976].)  Although  a  full  three-dimensional  magnetohydrodynamic 
treatment  of  this  expansion  process  is  beyond  the  scope  of  the  present  paper 
(see  § I V  for  further  discussion),  to  order  of  magnitude  one  may  assume  that  a 

disturbance  of  the  surrounding  field  lines  proceeds  horizontally  at  the  ion 

1/  - 

sound  speed  cg  »  (kTe/mH)  2  ,  which  is  approximately  equal  to  the  Alfven  speed 

»  B/(4»nmH)//2  since  the  gas  and  magnetic  pressures  are  comparable  (Fig¬ 
ure  1).  After  a  time  T  the  disturbance  encounters  the  neighboring  loop  Lj 
(Figure  2),  allowing  L2  to  release  Its  stored  free  energy  (in  the  form  of 


1 


1 


currents)  through  interaction  with  the  disturbed  field  lines  in  a 
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neutral-sheet-type  configuration  (see  Sturrock  1968),  and  giving  rise  to 
another  burst  of  hard  X-ray  emission,  spatially  distinct  from  the  location  of 
the  first  burst  (see  §1). 

In  order  for  the  phenomenon  just  described  to  explain  the  features  of 
EFB's,  the  interaction  time  T  must  be  roughly  equal  to  the  temporal  separation 
tEFB  EFB'S*  Thus,  the  lateral  separation  D  of  the  loops  L1  and  L2  must  be 
of  order 


D 


V  T 
A  EFB 


(28) 


which,  setting  t  »  5  -  10  s  (de  Jager  and  de  Jonge  1978)  and  VB  "5  x  10^ 
EFB  A 

cm  s'1  (using  B  “  500  G  and  the  density  values  shown  in  Figure  1),  gives 


D 


(3  -  5)  x  108 


cm  , 


(29) 


an  entirely  reasonable  value. 

The  situation  envisaged  ("driven  merging  flux")  is  somewhat  similar  to 
that  in  the  flare  model  by  Gold  and  Hoyle  (1960),  or  in  the  emerging  flux 
model  of  Heyvaerts,  Priest,  and  Rust  (1977).  However,  there  are  some  signif¬ 
icant  differences,  as  Figure  2  shows.  Because  the  two  loops  under  considera¬ 
tion  are  part  of  the  same  active  region  complex,  which  we  assume  to  have  a 
simple  bipolar  structure,  the  polarities  of  the  corresponding  footpoints  of 
the  two  loops  are  the  same ;  such  a  situation  gives  no  reversal  of  the  toroidal 
(i.e.  longitudinal)  component  of  the  magnetic  field  at  the  contact  plane  be- 

I 

tween  and  L2  and  so  this  toroidal  component  of  the  field  cannot  be  respon¬ 
sible  for  any  subsequent  energy  release  as  in  the  model  of  Heyvaerts,  Priest, 


r 


i 
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and  Rust  (1977).  Further,  the  (poloidal)  currents  corresponding  to  these  tor¬ 
oidal  fields  are  antiparallel  in  the  contact  plane  and  therefore  repel,  in¬ 
stead  of  attracting  in  the  (field-reversed)  Gold  and  Hoyle  (1960)  model.  How¬ 
ever,  by  the  same  considerations,  any  twist  (i.e.,  departure  from  a  potential 
configuration)  of  the  magnetic  field  lines  will  be  in  the  same  sense  as  viewed 
by  a  remote  observer,  which  leads  to  a  reversal  of  the  poloidal  components  of 
the  magnetic  field  in  the  contact  plane,  providing  a  situation  suitable  for 
reconnection  of  this  poloidal  component  in  a  neutral-sheet-type  scenario  (see 
Sturrock  1968).  The  force  required  to  overcome  the  repulsion  of  the  poloidal 
currents  is  the  hydrodynamic  pressure  of  the  heated  transition  region  in  the 
flaring  loop;  a  simple  calculation  shows  that  the  condition  B  “  1  is  equiva¬ 
lent  to  the  condition  that  the  gas  pressure  gradient  can  overcome  the  electri¬ 
cal  repulsion  of  these  poloidal  currents. 

This  "poloidal  flux  annihilation"  aspect  of  the  present  model  gives,  upon 
further  consideration,  some  rather  pleasing  results.  First,  since  the  toroi¬ 
dal  component  of  the  magnetic  field  remains  intact  after  the  reconnection  pro¬ 
cess,  the  basic  spatial  configuration  of  the  loops  remain  intact,  permitting 
the  formation  of  the  soft  X-ray  emitting  post-flare  loops  frequently  observed 
(e.g.,  Pallavicini,  Serio,  and  Vaiana  1977).  Second,  during  the  reconnection 
3B/3t,  and  so  7  x  E,  is  in  the  poloidal  direction.  This  gives  an  induced  E 
vector  which  is  predominantly  toroidal,  i.e.,  parallel  to  the  remaining  mag¬ 
netic  field  lines,  thus  allowing  efficient  acceleration  of  particles  along 


! 


these  field  lines.  (In  standard  neutral-sheet-type  models  [e.g.,  Sturrock 
1968]  3B/3t  and  7  x  E  are  mainly  vertical,  resulting  in  an  induced  E  vector 
which  is  horizontal,  i.e.,  perpendicular  to  the  remaining  [closed]  field 
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lines.)  Third,  there  is  no  need  for  a  complex  magnetic  topology  at  the  loop 

footpoints,  such  as  is  required  by  both  the  Gold  and  Hoyle  (1960)  model  and 

the  emerging  flux  model  of  Heyvaerts,  Priest,  and  Rust  (1977). 

One  may  also  estimate  the  energy  released  in  the  annihilation  of  this 

poloidal  flux.  Assuming  a  twist  of  *  if  in  each  flux  tube  (see  Barnes  and 

Sturrock  1972),  a  toroidal  field  strength  of  some  300  G  (see  above),  and  a 
2  17  2 

footpoint  area  A  *  D  “3x10  cm  ,we  find  that  the  total  energy  released 
is 


2  2  2 

e  B  D  L  fD)2 
Erel  "  8r  ' 


(30) 


where  L  is  the  length  of  the  flux  tube  (“  3  x  109  cm)  and  e  is  the  fraction  of 
the  poloidal  field  that  is  destroyed.  Now,  the  proposed  model  relies  on  the 
fact  that  e  <  1,  since  some  poloidal  flux  must  remain  after  each  loop  "encoun¬ 
ter"  if  the  "domino"  effect  leading  to  successive  EFB’s  is  to  work.  Thus, 
setting  £  “  1/3,  B  =  500  G,  D  =  5  x  10^cm,  and  L  =  3  x  10®  cm,  we  obtain 


28 

E  ,  -  3  x  10  ergs  .  (31) 

rel 

We  may  compare  this  with  the  energy  required  to  produce  a  single  hard  X-ray 
EFB.  Hard  X-ray  spectra  in  moderately  large  flares  may  be  represented  well  by 
the  power-law  form 


-v  -2  -1  -1 

I(£)  *  a£  photons  cm  s  keV  , 


(32) 


with  a  “  10  and  y  (“  6  -1)  “3-4  (see  Hoyng,  Brown,  and  van  Beek  1976).  To 
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produce  such  a  burst  by  thick  target  bremsstrahlung  of  a  beam  of  high  energy 
electrons  requires  an  injection  rate  of  electrons  with  energies  >  25  keV  of 

F25  =  4  x  1033a(Y-1)2B(Y->/2  ,V2  )(25)"Y  *  1035-1036electrons  s"1  (33) 

(Hoyng,  Brown,  and  van  Beek  1976),  corresponding  to  a  total  injected  power 

Y  27  2R  —1 

p  =  F25  X  25  keV  =  3  x  10  -  3  x  10  erg  s  .  (34) 

Comparing  eqs.  (31)  and  (34)  with  a  FWHM  for  the  EFB  of  “5s  (de  Jager  and  de 
Jonge  1978)  indicates  agreement  to  within  an  order  of  magnitude,  which  consid 
ering  the  crudeness  of  the  above  argument,  is  quite  acceptable. 

This  completes  the  scenario  of  the  interacting  loop  model,  whose  energet 
ics  and  timescales  have  been  shown  to  agree  well  with  observations.  In  the 
final  section  we  shall  discuss  some  features  of  the  model  in  greater  detail, 
and  also  point  out  areas  for  future  work. 

IV.  DISCUSSION 

The  analysis  of  the  preceding  two  sections  is  of  course  somewhat  schema¬ 
tic,  and  is  intended  only  to  provide  a  general  description  of  the  main  fea¬ 
tures  of  the  model.  Nevertheless,  it  not  only  accounts  for  the  EFB  descrip¬ 
tion  of  hard  X-ray  bursts,  but  also,  for  plausible  parameters,  successfully 
accounts  for  both  the  observed  intensity  and  temporal  spacing  of  EFB ' s.  The 
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only  a  priori  requirements  of  the  model  are  (i)  the  initial  energy  release  in 
one  of  the  active  region  loops  (L^  in  Figure  2),  which  may  occur  by  any  of  the 
currently  proposed  mechanisms  (e.g.,  the  tearing  mode  instability  in  a 
stressed  arch  geometry—Spicer  1976,  1977;  van  Hoven  1979),  and  (ii)  the  exis¬ 
tence  of  a  large  scale  current  pattern  in  the  active  region,  to  provide  the 
poloidal  flux  component  which  is  annihilated  in  the  "domino”  process  of  suc¬ 
cessive  energy  releases. 

Clearly  many  aspects  of  the  model  require  more  detailed  study.  For  exam¬ 
ple,  in  Figure  2  it  is  tacitly  assumed  that  the  energized  flux  tube  expands 
along  its  entire  length  as  a  result  of  the  gas  pressure  forces  which  are 
localized  at  its  base.  Our  intuitive  justification  for  this  is  that  the  ten¬ 
sion  in  the  magnetic  lines  of  force  will  tend  to  "straighten  out"  any  "bulge” 
formed  by  a  local  expansion  of  the  field;  however,  the  timescale  for  the  en¬ 
tire  loop  to  respond  to  this  tensile  force  will  be  of  order  L/VA  corona,  which 
for  very  large  (L/D)  may  be  larger  than  the  inter-loop  collision  timescale  tc 
(eq.  [29]),  even  allowing  for  the  reduced  density,  and  corresponding  higher 
Alfven  velocity,  in  the  corona.  A  full  three-dimensional  magnetohydrodynamic 
analysis  of  the  situation  is  clearly  necessary  to  resolve  this  issue  satisfac¬ 
torily,  just  as  it  is  also  required  to  adequately  model  the  expansion  of  the 
flux  tube  in  the  (difficult)  0  “  1  regime  (see  remarks  at  the  beginning  of 
§111). 

In  addition  to  these  issues,  the  details  of  the  interaction  between  "col¬ 
liding"  flux  tubes  (see  Figure  2)  needs  to  be  explored  more  fully.  This  re¬ 
quires  study  of  driven  magnetic  reconnection  processes,  such  as  have  been  dis¬ 
cussed  previously,  for  the  case  of  flux  driven  by  magnetic  buoyancy,  by 
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Heyvaerts,  Priest,  and  Rust  (1977),  and  by  references  therein.  This  will  en¬ 
able  us  to  better  determine  the  temporal  and  spatial  characteristics  of  the 
energy  release  in  the  quiescent  loop  (loop  Lj  in  Figure  2). 

I  thank  J.  Leach  for  stimulating  discussions  which  led  to  the  ideas  de¬ 
veloped  in  this  paper,  A.  N.  McClymont  for  useful  discussions  regarding  the 
hydrodynamics  of  the  solar  atmosphere  during  flares,  and  S.  K.  Antiochos  and 
P.  A.  Sturrock  for  their  comments  on  the  manuscript.  The  author  has  benefited 
from  discussions  held  at  the  Study  of  Energy  Release  in  Flares  (SERF)  Workshop 
held  at  Stanford  in  August,  1980?  this  workshop  was  sponsored  by  the  National 
Science  Foundation  and  by  the  Scientific  Committee  on  Solar-Terrestrial 
Physics  (SCOSTEP).  Financial  support  from  grants  NASA  NGL  05-020-272,  NASA 
NAGW-92,  and  contract  ONR  NOOO 14-75-C-0673,  is  also  gratefully  acknowledged. 
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FIGURE  CAPTIONS 


Fig.  1.  -  Variation  of  density  n  (cm-3)  and  electron  temperature  T@  (K)  with 
-2 

column  density  N  (cm  ),  from  the  Somov,  Syrovatskii,  and  Spektor  (1981) 

calculations  of  the  hydrodynamic  response  of  the  solar  chromosphere  to  an 

energy  input  in  the  form  of  a  beam  of  non-thermal  electrons.  The  preflare 

values  of  n  and  Tg  are  shown  dashed;  see  §11  for  a  discussion  of  these  initial 

conditions.  The  values  in  the  flare  correspond  to  a  time  5  seconds  after  the 

start  of  the  heating  and  correspond  to  the  largest  enhancements  in  their 

calculations.  The  ion  temperature  T^  (not  shown)  is  much  less  than  Tg.  Also 

shown  is  the  pressure  P  =  nTg  (cm-3  K)  and  the  equipartition  magnetic  field 

B  (gauss;  see  eq.  [3] ).  Note  the  large  values  of  P  and  B  around  the 
“4  e4 

transition  zone,  caused  by  the  impulsive  heating  of  plasma;  the  value  of  Beg 
there  is  of  order  100  G. 


Fig.  2.  -  Schematic  of  the  flux  tube  interaction.  Loops  and  L2  are  typical 
members  of  the  same  bipolar  active  region  complex  and  lie  across  the  neutral 
line  (heavy  dashed)  as  shown.  Their  magnetic  fields  are  also  twisted  in  the 
same  sense  due  to  a  global  active  region  current  field  J.  At  some  time  t  =  0 
loop  L1  flares  (the  energy  release  being  signified  by  the  star  at  the  loop 
apex),  causing  energy  to  flow  downwards  toward  the  chromosphere  (the  energy 
transport  mechanism  is  here  taken  to  be  non-thermal  electron  bombardment,  so 
that  we  may  analyze  its  effects  in  some  detail  -  see  §11  -  however,  other 
transport  mechanisms  [see  §11]  will  have  the  same  qualitative  effect  on  the 


evolution  of  the  system).  Upon  interaction  with  the  chromosphere,  the  trans- 
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ported  energy  heats  the  plasma  to  very  high  (■  107  K)  temperatures,  causing  a 
large  pressure  enhancement  at  this  point  (see  Figure  1).  This  causes  loop  L. 
to  expand,  driving  a  horizontal  magnetohydrodynamic  disturbance  (small  hori¬ 
zontal  arrows  at  the  loop  base).  At  some  later  time  the  surface  of  the 
disturbed  field  lines  has  the  appearance  1  (right  figure;  the  original  loop 
L.j  is  shown  dashed),  and  is  in  contact  with  loop  Lj'  with  the  poloidal  com- 

I 

ponents  of  the  magnetic  fields  of  and  L2  antiparallel  in  the  contact 
plane.  These  components  reconnect  in  a  neutral-sheet-type  scenario  and  cause 
an  impulsive  energy  release  in  the  hitherto  quiescent  loop  L2» 

For  the  purpose  of  the  illustration,  the  poloidal  component  of  the  mag¬ 
netic  fields  has  been  exaggerated.  Also,  loop  '  has  been  constructed  as  a 
simple  radial  dilation  of  loop  I,.';  this  is  not  necessarily  justified  in  prac¬ 
tice  (see  §IV) . 
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